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What is Claimed is: 

1. A method executed in a computer system for producing a combined system of partial 
differential equations comprising: 

representing each of a plurality of systems as an application mode modeling physical 
quantities of said each system; 

determining a representation of a partial differential equation system for each appKcation 
mode corresponding to one of said plurality of systems using at least one non-local coupling, 
said at least one non-local coupling determining a value in at least one point depending on a 
value from at least one other point; and 

forming said combined system of partial differential equations using partial differential 
equation systems associated with said plurality of systems. 

2. The method of Claim 1 wherein said non-local variable defines a value from a first 
domain in a first geometry to another domain in a second geometry. 

3. The method ofClaim 2, fiirther comprising: 

forming, for each of said first and second geometries, a system of partial differential 
equations each having associated coupling variables. 

4. The method of Claim 2, wherein at least one of said partial differential equation 
systems uses at least one local coupling. 
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5. The method of Claim 2, wherein said first and second geometries are the same. 

6. The method of Claim 2, wherein said first and second geometries are different. 

7. The method of Claim 1, further comprising: 

defining a non-local coupling wherein a value of a quantity on a boundary of a first 
domain is referenced on parallel lines extending into said first domain. 

8. The method of Claim 2, fiirther comprising: 

defining a non-local coupling in which a boundary condition associated with said first 
domain is defined using a value of an integral over a portion of one of: said first domain and said 
second domain. 

9. The method of Claim 1, fiarther comprising: 

defining a non-local coupling using at least one of: a mapped variable and an integrated 
variable. 

10. The method of Claim 4, fiirther comprising: 

defining a local coupling using at least one of: a basic variable, an auxiliary variable, and 
a glued variable. 
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11. The method of Claim 1 , further comprising: 

defining a non-local coupling variable using at least one of: an extrusion variable, a 
projection variable, and a scalar coupling variable. 

12. The method of Claim 1, fiulher comprising: 

determining a stiffness matrix by determining at least one of a Jacobian of a variable and 
a value of a variable in accordance with a type of said variable wherein said combined system of 
partial differential equations is in weak form, said stiffiiess matrix being a Jacobian matrix of a 
residual vector with respect to a number of degrees of freedom, said Jacobian of a variable being 
represented as at least one contribution determined in accordance with a number of degrees of 
freedom; and 

determining said residual vector by determining at least one of a Jacobian of a variable 
and a value of a variable in accordance with a type of said variable wherein said combined 
system of partial differential equations is in weak foim. 

13. The method of Claim 12, fiirther comprising: 

converting said combined system of partial differential equations from general form to 
weak form. 
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14. The method of Claim 12, wherein said detemiining said stiffhess matrix further 
comprises: 

determining at least one of a Jacobian of a variable and a value of a variable in 
accordance with points included in a quadrature formula and with other points in accordance 
with coupUng variables. 

15. The method of Claim 14, wherein said determining said residual vector further 
comprises: 

determinmg at least one of a Jacobian of a variable and a value of a variable in 
accordance with pomts included in a quadrature formula and with other points in accordance 
with coupling variables. 

16. The method of Claim 1, further comprising: 

determining a value of a variable in accordance with a type of said variable used in at 
least one of said partial differential equation systems. 

17. The method of Claim 16, wherein variables are recursively evaluated in accordance 
with variable type. 

1 8. The method of Claim 17, wherein said determining a value of a variable in 
accordance with a type is used in determining at least one of a: stiffiiess matrix, constraint 
matrix, residual vector and a constraint residual vector. 
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19. The method of Claim 18, wherein said type is one of: a basic variable, an auxihary 
variable, a glued variable, a mapped variable, and an integrated variable. 

20. The method of Claim 1, further comprising: 

determining a Jacobian of a variable in accordance with a type of said variable used in at 
least one of said partial differential equations wherein said Jacobian of a variable is represented 
in accordance with a number of degrees of freedom. 

21. The method of Claim 20, wherein a Jacobian of a variable is recursively evaluated in 
accordance with variable type. 

22. The method of Claim 2 1 , wherein said determining a Jacobian of a variable in 
accordance with a type is used in determining at least one of: a stiffeess matrix, a residual vector, 
constraint residual vector, and a constraint matrix. 

23. The method of Claim 22, wherein said type is one of: a basic variable, an auxihary 
variable, a glued variable, a mapped variable, and an integrated variable. 

24. The method of Claim 1 , wherein said at least one non-local coupling includes a 
variable having a dependency on another variable at at least one distant point. 

25. The method of Claim 24, wherein said other variable is in the same geometry as said 
variable. 
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26. The method of Claim 24, wherein said other variable is in a different geometry from 
said variable. 

27. The method of Claim 4, wherein said local coupKng includes a variable having a 
dependency only on values of other variables at the same point. 

28 . The method of Claim 1 , further comprising: 

defining a non-local coupling used in at least one of a: subdomain, boundary, edge, and 
point that obtains a value at one of: a subdomain, boundary, edge, and point. 

29. The method of Claim 4, further comprising: 

defining a local coupling using at least one of: an expression variable and a boundary 
coupled variable. 

30. The method of Claim 1, fiirther comprising: 

defining a non-local coupling wherein a value of an integral of a variable along parallel 
Hnes extending into a furst domain is used on a boundary on said first domain. 



31. The method of Claim 1, further comprising: 
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defining a boundary condition on one boundary in terms of a value of a variable on 
another boundary wherein said value is mapped in accordance with a coordinate transformation. 

32. The method of Claim 1, further comprising: 

defining a boundary condition in terms of a variable defined at a single point. 
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33. A method executed in a computer system for assembling a finite element 
discretization of a system of weak partial differential equations comprising: 

determining a stif&ess matrix by evaluating at least one of a Jacobian of a variable and a 
value of a variable in accordance with a type of said variable included in said system, said 
Jacobian of said variable being represented as at least one contribution in accordance with a 
number of degrees of freedom; and 

determining a residual vector by evaluating at least one of a Jacobian of a variable and a 
value of a variable in accordance with a type of said variable included in said system, said 
Jacobian of said variable being represented as at least one contribution determined in accordance 
with a number of degrees of freedom. 

34. The method of Claim 33, wherein said system includes at least one non-local 
coupling and said method fiirther includes evaluating at least one variable in said non-local 
coupling in at least one of: said determining said stiffness matrix and said determining said 



IS residua] vector. 



35. The method of Claim 33, wherein said system includes at least one local coupling 
and said method further includes evaluating at least one variable in said local coupling in at least 
one of: said determining said stiffiiess matrix and said determining said residual vector. 
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36. The method of Claim 33, further comprising: 

determining a constraint matrix by evaluating at least one of said Jacobiati of a variable 
and a value of a variable in accordance with a type of said variable included in said system 
wherein said Jacobian of said variable is represented as at least one contribution in accordance 
with a number of degrees of freedom. 

37. The method of Claim 36, further comprismg: 

determining a constraint residual vector by determining a value of at least one variable 
included in said system in accordance with a type of said at least one variable. 

38. The method of Clahn 37, wherein said type is one of: a basic variable, an auxiUary 
vaiable, a glued variable, a mapped variable, an integrated variable. 

39. The method of Claim 33, wherein said type is one of: a basic variable, an auxihary 
variable, a glued variable, a mapped variable, an integrated variable. 
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40. A computer program product for producing a combined system of partial differential 
equations comprising machine executable code for: 

representing each of a plurality of systems as an apphcation mode modehng physical 
quantities of said each system; 

determining a representation of a partial differential equation system for each application 
mode corresponding to one of said plurality of systems using at least one non-local couphng, 
said at least one non-local coupling determining a value in at least one point depending on a 
value from at least one other point; and 

forming said combined system of partial differential equations using partial differential 
equation systems associated with said plurality of systems. 

41 . The computer program product of Claim 40, wherein said non-local variable defines a 
value from a first domain in a first geometry to another domain in a second geometry. 

42. The computer program product of Claim 41 , further comprising machine executable 
code for: 

forming, for each of said first and second geometries, a system of partial differential 
equations each having associated coupling variables. 

43. The computer program product of Claim 41, wherein at least one of said partial 
differential equation systems uses at least one local coupling. 
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44. The computer program product of Claim 41, wherein said first and second 
geometries are the same. 

45. The computer program product of Claim 41, wherein said first and second 
geometries are different. 

46. The computer program product of Claim 40, fiirther comprising machine executable 
code for: 

defining a non-local coupling wherein a value of a quantity on a boundary of a first 
domain is referenced on parallel Imes extending into said first domain. 

47. The computer program product of Claim 41, fiuther comprising machine executable 
code for: 

defining a non-local coupling in which a boundary condition associated with said first 
domain is defined using a value of an integral over a portion of one of: said first domam and said 
second domain. 

48. The computer program product of Claim 40, fiirther comprising machine executable 
code for: 

defining a non-local coupling using at least one of: a mapped variable and an integrated 
variable. 
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49. The computer program product of Claim 43, further comprising machine executable 
code for: 

defining a local coupling using at least one of: a basic variable, an auxiliary variable, and 
a glued variable. 

50. The computer program product of Claim 40, further comprising machine executable 
code for: 

defining a non-local coupling variable using at least one of: an extrusion variable, a 
projection variable, and a scalar coupling variable. 

5 1 . The computer program product of Claim 40, fiuther comprising machine executable 
code for: 

determining a stif&iess matrix by determining at least one of a Jacobian of a variable and 
a value of a variable in accordance with a type of said variable wherein said combined system of 
partial differential equations is in weak form, said stiffiiess matrix being a Jacobian matrix of a 
residual vector with respect to a number of degrees of freedom, said Jacobian of a variable being 
represented as at least one contribution determined in accordance with a number of degrees of 
freedom; and 

determining said residual vector by determining at least one of a Jacobian of a variable 
and a value of a variable in accordance with a type of said variable wherein said combined 
system of partial differential equations is in weak form. 
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52. The computer program product of Claim 51, further comprising machine executable 
code for: 

converting said combined system of partial differential equations from general form to 
weak form. 

53. The computer program product of Claim 51, wherein said machine executable code 
for determining said stiffiiess matrix further comprises: 

machine executable code for determining at least one of a Jacobian of a variable and a 
value of a variable in accordance with points included in a quadrature formula and with other 
points in accordance with coupling variables. 

54. The computer program product of Claim 53, wherein said machine executable code 
for determining said residual vector further comprises: 

machine executable code for determining at least one of a Jacobian of a variable and a 
value of a variable in accordance with points included in a quadrature formula and with other 
points in accordance with coupling variables. 

55. The computer program product of Claim 40, further comprising machine executable 
code for: 

determining a value of a variable in accordance with a type of said variable used in at 
least one of said partial differential equation systems. 
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56. The computer program product of Claim 55, wherein variables are recursively 
evaluated in accordance with variable type. 

57. The computer program product of Claim 56, wherein said machine executable code 
for determining a value of a variable in accordance with a type is used in machine executable 
code for determining at least one of a: stif&ess matrix, constraint matrix, residual vector and a 
constraint residual vector. 

58. The computer program product of Claim 57, wherein said type is one of: a basic 
variable, an auxiliary variable, a glued variable, a mapped variable, and an integrated variable. 

59. The computer program product of Claim 40, further comprising machine executable 
code for: 

determining a Jacobian of a variable in accordance with a type of said variable used in at 
least one of said partial differential equations wherein said Jacobian of a variable is represented 
in accordance with a number of degrees of freedom. 

60. The computer program product of Claim 59, wherein a Jacobian of a variable is 
recursively evaluated in accordance with variable type. 
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61 . The computer program product of Claim 60, wherein said machine executable code 
for determining a Jacobian of a VMiable in accordance with a type is used in machine executable 
code for determining at least one of: a stiffness matrix, a residual vector, constraint residual 
vector, and a constraint matrix. 

62. The computer program product of Claim 6 1 , wherein said type is one of: a basic 
variable, an auxiUary variable, a glued variable, a mapped variable, and an integrated variable. 

63. The computer program product of Claim 40, wherein said at least one non-local 
coupling includes a variable having a dependency on another variable at at least one distant 
point. 

64. The computer program product of Claim 63, wherein said other variable is in the 
same geometry as said variable. 

65. The computer program product of Claim 63, wherein said other variable is in a 
different geometry from said variable. 

66. The computer program product of Claim 43, wherein said local coupling includes a 
variable having a dependency only on values of other variables at the same point. 
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67. The computer program product of Claim 40, further comprising machine executable 
code for: 

defining a non-local coupling used in at least one of a: subdomain, boundary, edge, and 
point that obtains a value at one of: a subdomain, boundary, edge, and point. 

68. The computer program product of Claim 43, further comprising machine executable 
code for: 

defining a local coupling using at least one of: an expression variable and a boundary 
coupled variable. 

69. The computer program product of Claim 40, further comprising machine executable 
code for: 

defining a non-local coupling wherem a value of an mtegral of a variable along parallel 
lines extending into a first domain is used on a boundary on said first domain. 

70. The computer program product of Claim 40, fiirther comprising machine executable 
code for: 

defining a boundary condition on one boundary in terms of a value of a variable on 
another boundary wherein said value is mapped m accordance with a coordinate transformation. 
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71. The computer program product of Claim 40, further comprising machine executable 
code for: 

defining a boundary condition in terms of a variable defined at a single point. 
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72. A computer program product for assembling a finite element discretization of a 
system of weak partial differential equations comprising machine executable code for: 

determining a stiffness matrix by evaluating at least one of a Jacobian of a variable and a 
value of a variable in accordance with a type of said variable included in said system, said 
Jacobian of said variable being represented as at least one contribution in accordance with a 
number of degrees of freedom; and 

determining a residual vector by evaluating at least one of a Jacobian of a variable and a 
value of a variable in accordance with a type of said variable included in said system, said 
Jacobian of said variable being represented as at least one contribution determined in accordance 
with a number of degrees of freedom. 

73. The computer program product of Claim 72, wherein said system includes at least 
one non-local coupling and said computer program product further includes machine executable 
code for evaluating at least one variable in said non-local coupling in at least one of: said 
machine executable code for determining said stiffness matrix and said machine executable code 
for determining said residual vector. 

74. The computer program product of Claim 72, wherein said system includes at least 
one local coupling and said computer program product further includes machine executable code 
for evaluating at least one variable in said local coupling in at least one of: said machine 
executable code for determining said stiffiiess matrix and said machine executable code for 
determining said residual vector. 
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75. The computer program product of Claim 72, further comprising machine executable 
code for: 

determining a constraint matrix by evaluating at least one of said Jacobian of a variable 
and a value of a variable in accordance with a type of said variable included in said system 
wherein said Jacobian of said variable is represented as at least one contribution in accordance 
with a number of degrees of freedom. 

76. The computer program product of Claim 75, further comprising machine executable 
code for: 

determining a constraint residual vector by determining a value of at least one variable 
included in said system in accordance with a type of said at least one variable. 

77. The computer program product of Claim 76, wherein said type is one of: a basic 
variable, an auxiliary variable, a glued variable, a mapped variable, an integrated variable. 

78. The computer program product of Claim 72, wherein said type is one of: a basic 
variable, an auxiUary variable, a glued variable, a mapped variable, an mtegrated variable. 
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